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ABSTRACT 

A  laser-doppler  anemometer  was  used  to  obtain  flow-field 
velocity  measurements  in  a  4:1  pressure  ratio,  4.54  kg/s  (10  Ibm/s), 
centrifugal  impeller,  with  splitter  blades  and  backsweep,  which  was 
configured  with  a  vaneless  diffuser.  Measured  through-flow  velocities 
are  reported  for  ten  quasi-orthogonal  survey  planes  at  locations  ranging 
from  1%  to  99%  of  main  blade  chord.  Measured  through-flow 
velocities  are  compared  to  those  predicted  by  a  3-D  viscous  steady 
flow  analysis  (Dawes)  code. 

The  measurements  show  the  development  and  progression  through 
the  impeller  and  vaneless  diffuser  of  a  through-flow  velocity  deficit 
which  results  from  the  tip  clearance  flow  and  accumulation  of  low 
momentum  fluid  centrifuged  from  the  blade  and  hub  surfaces.  Flow 
traces  from  the  CFD  analysis  show  the  origin  of  this  deficit  which 
begins  to  grow  in  the  inlet  region  of  the  impeller  where  it  is  first 
detected  near  the  suction  surface  side  of  the  passage.  It  then  moves 
toward  the  pressure  side  of  the  channel,  due  to  the  movement  of  tip 
clearance  flow  across  the  impeller  passage,  where  it  is  cut  by  the 
splitter  blade  leading  edge.  As  blade  loading  increases  toward  the  rear 
of  the  channel  the  deficit  region  is  driven  back  toward  the  suction 
surface  by  the  cross-passage  pressure  gradient.  There  is  no  evidence 
of  a  large  wake  region  that  might  result  from  flow  separation  and  the 
impeller  efficiency  is  relatively  high.  The  flow  field  in  this  impeller  is 
quite  similar  to  that  documented  previously  by  NASA  Lewis  in  a  large 
low-speed  backswept  impeller. 

NOMENCLATURE 

N  =  rotational  speed,  rev/min 

^design  =  design  rotational  speed,  rev/min 

Ncpp  =  number  of  measurement  bins  per  passage 


PS  =  pressure  surface  of  a  blade 
SS  =  suction  surface  of  a  blade 
Utip  =  tip  velocity  at  impeller  exit,  m/s 
Vgbs  =  absolute  velocity  in  the  laser  measurement  plane,  m/s 
Vqn,  =  velocity  component  normal  to  beam  path  in  the  meridional 
plane.  Also  referred  to  as  through-flow  velocity,  m/s 
Ve  =  velocity  component  in  direction  of  rotation,  m/s 
x'  =  through-flow  direction  in  the  laser  measurement  plane 
6  =  circumferential  direction 

P  =  absolute  flow  angle  in  measurement  plane,  P=  tan'^(Ve  /VqJ 
pjjjjj,  =  minimum  flow  angle  at  each  downstream  radius  ratio  and  span 

INTRODUCTION 

Laser  anemometry  has  been  used  by  several  investigators  to 
obtain  measurements  of  the  flow  fields  within  high  speed  centrifugal 
compressors  over  the  last  20  years.  The  first  reported  measurements 
were  obtained  by  Eckardt  (1976)  using  a  Laser-2-Focus  Velocimeter. 
Since  then,  several  other  researchers  (  Adler  and  Levy,  1979;  Eckardt, 
1979;  Krain,  1981,  1987;  Ahmed  and  Elder,  1990;  Pak  et  al.,  1993; 
Krain  et  al.,  1995)  have  used  laser  anemometers  to  measure  the  flow 
fields  within  high-speed  centrifugal  impellers  configured  with  both 
vaned  and  vaneless  diffusers.  Hathaway  et  al.  (1993)  reported  results 
from  a  study  in  a  large  low-speed  compressor. 

In  the  previous  studies,  a  three-dimensional  flow  field  has  been 
observed  in  the  impeller.  The  measured  flows  have  been  characterized 
by  a  region  of  fluid  with  low  through-flow  velocity  which  is  first 
observed  at  various  meridional  and  pitchwise  locations.  This  region  of 
low  velocity  fluid  has  generally  been  attributed  to  a  separation 
initiated  by  blade  and  shroud  curvature.  At  the  impeller  discharge,  the 
flow  profile  ranges  from  a  jet-wake  type  flow  to  a  fairly  uniform  flow 
as  reported  by  Adler  and  Levy  (1979).  The  radially  ending  impellers 
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had  a  discharge  profile  resembling  the  jet-wake,  whereas  the  backswept 
impellers  demonstrated  a  more  uniform  profile.  The  depth  and  location 
of  the  wake  zone  appears  to  have  been  modified  by  the  amount  of 
blade  backsweep  that  was  applied  at  the  impeller  exit. 

Hathaway  et  al.  (1993)  obtained  detailed  measurements  of  the 
secondary  flow  fields  in  a  large-low-speed  impeller  and  compared 
them  to  the  results  of  a  3D  Navier-Stokes  simulation  using  the 
BTOB3D  code  (Dawes,  1988).  This  work  suggested  that  the  origin  of 
the  low  momentum  fluid  may  lie  in  a  vortex  which  is  formed  by  the 
interaction  of  a  tip  clearance  jet  with  the  centrifuged  suction  surface 
boundary  layer.  The  vortex  first  migrates  toward  the  pressure  surface 
and  then  toward  the  suction  surface.  Krain  (1987)  showed  that  the 
downstream  regions  of  low  momentum  fluid  have  the  characteristics 
of  a  vortex  flow. 

The  purpose  of  the  work  reported  herein  was  two-fold.  First,  to 
obtain  a  detailed  survey  of  the  flow  fields  in  a  high  speed  centrifugal 
impeller  with  splitters  so  as  to  obtain  a  data  base  for  validation  of  the 
CFD  codes  that  can  analyze  such  geometries.  Second,  to  see  if 
evidence  of  the  vortex  initiation  identified  in  the  large  low-speed 
compressor  could  also  be  found  in  a  smaller  high-speed  machine. 

The  impeller  used  in  this  study  was  configured  with  both 
vaneless  and  vaned  diffusers.  The  overall  performance  of  each 
configuration  is  presented.  Laser  surveys  were  completed  in  each 
configuration,  but  only  those  obtained  in  the  impeller  when  it  was 
configured  with  the  vaneless  diffuser  are  presented  here.  To  validate 
the  laser  measurements,  velocity  profiles  calculated  from  total  pressure 
measurements  in  the  vaneless  diffuser  are  presented  and  compared  to 
those  measured  using  laser  anemometry.  The  measurements  are 
supported  further  by  a  BTOB3D  code  (Dawes,  1988)  analysis  which 
predicts  a  flow  field  that  is  similar  to  the  one  that  was  measured  and 
shows  the  origin  of  the  flow  features  that  were  observed. 

TEST  COMPRESSOR 

The  test  compressor  is  a  Allison  Engine  Company  design  that  was 
scaled  up  to  a  flow  size  of  4.54  kg/s  (10  Ibm/s)  from  the  original  size 
of  1.66  kg/s  (3.655  Ibm/s).  Complete  aerodynamic  and  mechanical 
design  descriptions,  including  impeller  and  diffuser  geometries,  are 
given  by  McKain  and  Holbrook  (1982). 

The  stage  (impeller  with  vaned  diffuser)  was  designed  to  produce 
a  pressure  ratio  of  4:1  at  the  design  mass  flow.  The  standard  day 
corrected  speed  for  the  design  flow  condition  is  21,789  rev/min  with 
an  exit  tip  speed  of  492  m/s  (1615  ft/s).  The  inlet  relative  Mach 
numbers  on  the  suction  surface  range  from  0.9  at  the  tip  to  0.45  at  the 
hub.  The  impeller  contains  15  main  blades  with  15  splitter  blades  and 
has  50  degrees  of  backsweep  from  radial  at  the  discharge.  The  splitter 
blade  leading  edge,  located  at  30  percent  of  main  blade  chord,  is  offset 
shghtly  toward  the  main  blade  suction  surface  in  order  to  produce  an 
even  flow  split.  Both  the  main  blades  and  splitter  blades  are  formed 
from  quasi-normal  straight-line  elements  between  the  hub  and  tip.  The 
inlet  tip  diameter  is  210  mm  (8.264  in)  and  the  inlet  blade  height  is 
64  mm  (2.501  in).  The  exit  diameter  is  431  mm  (16.986  in)  and  the 
exit  blade  height  is  17  mm  (0.671  in).  All  dimensions  are  for  hot 
conditions  at  100  percent  of  design  speed. 


Figure  1 :  Centrifugal  Compressor  Flow  Path  With  Laser 
Windows. 

A  test-rig  cross  section  showing  the  flow  path  through  the 
impeller  and  vaneless  diffuser  can  be  seen  in  Figure  1.  For  the 
impeller-with-vaneless-diffuser  configuration,  pressure  and  temperature 
rakes  were  located  at  measurement  station  A  which  was  at  a  radius 
ratio  of  1.18.  When  the  vaned  diffuser  was  installed  (stage 
configuration)  rakes  at  measurement  station  B  were  used. 

The  five  laser  windows  shown  in  Figure  1  are  from  one  of  four 
interchangeable  window  sets.  By  using  various  combinations  of 
windows  from  each  set,  it  was  possible  to  obtain  nearly  full  optical 
coverage  of  the  flow  path  through  the  impeller  and  vaneless  diffuser. 

The  test  facility  used  for  both  performance  measurements  and 
laser  anemometer  surveys  is  an  open  loop  facility  with  variable  inlet 
and  exhaust  pressure  control.  It  can  accommodate  both  axial  and 
centrifugal  compressors  in  a  wide  range  of  flow  sizes  and  pressure 
ratios.  A  complete  description  of  the  facility  is  given  by  Brokopp  and 
Gronski  (1992). 

COMPRESSOR  PERFORMANCE 

Aerodynamic  performance  maps  of  the  stage  and  impeller-with- 
vaneless-diffuser  configurations  are  shown  in  Figures  2  and  3, 
respectively.  Both  performance  maps  were  obtained  while  the  impeller 
was  running  with  a  tip  clearance  distribution  of  0.1524  mm  (0.006 
inch)  near  die  leading  edge,  0.61  mm  (0.024  inch)  near 
mid-chord,  and  0.203  mm  (0.008  inch)  near  the  exit. 

Stage  performance  is  based  upon  measurements  of  total 
temperature  and  total  pressure  that  were  obtained  from  fourteen  three 
and  four  element  rakes  spaced  evenly  at  measurement  station  B  (see 
Fig.  1).  These  were  compared  to  measurements  taken  upstream  of  the 
impeller  inlet  in  the  plenum  tank.  The  flow  rate  was  varied  between 
choke  and  near  surge  at  each  speed  using  a  downstream  throttle  valve. 
The  diffuser  controlled  the  choke  and  surge  mass  flows  at  all  speeds. 
The  lowest  flow  rate  shown  on  each  speed  line  in  Figure  2  is  the  last 
stable  operating  point  that  could  be  measured. 
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Figure  2:  Performance  Map  Of  4:1  Pressure  Ratio  Stage 


At  design  speed,  a  stage  total-total  peak  efficiency  of  0.846  was 
measured  at  a  mass  flow  of  4.39  kg/s  (9.70  Ibm/sec).  At  the  design 
flow  rate  and  speed,  a  total-total  efficiency  of  0.832  was  obtained  for 
a  total  pressure  ratio  of  3.97.  Variable  air  properties  were  used  in  the 
efficiency  calculation.  The  Mach  number  at  the  measurement  plane 
was  low  at  0.14,  so  the  measured  efficiency  compares  favorably  to  the 
0.830  total-static  efficiency  predicted  by  McKain  and  Holbrook  (1992) 
for  the  scaled  stage. 

Impeller-with-vaneless-diffuser  performance  is  shown  in  Figure  3. 
The  design  operating  point,  which  is  also  the  operating  point  used  for 
laser  surveys,  is  highlighted  in  the  figure.  Total  temperature  and 
pressure  were  determined  from  ten  evenly  spaced  three  and  four 
element  rakes  located  in  the  vaneless  diffuser  at  a  radius  ratio  of  1.18 
(location  A  in  Figure  1).  The  flow  range  at  each  speed  was  greater 
due  to  the  absence  of  the  vaned  diffuser.  Because  of  the  closer 
proximity  of  the  measurement  plane  to  the  impeller,  higher  pressures 
and  efficiencies  were  measured  at  all  speeds. 

At  design  speed,  a  peak  efficiency  of  0.867  occurred  at  the  design 
mass  flow  rate.  As  speed  was  reduced  to  70  percent  of  design,  peak 
efficiency  increased  to  0.910.  At  60  percent  speed  the  peak  efficiency 
decreased  to  0.897. 


Corrected  Mass  Flow  (Ibm/sec) 

Figure  3:  Performance  Map  Of  4:1  Pressure  Ratio  Impeller 
With  Vaneless  Diffuser 


LASER  ANEMOMETER  MEASUREMENT  TECHNIQUE 
General  Description 

A  single  component  laser-doppler  anemometer  operating  in  the 
back-scatter  mode  without  frequency  shifting  was  used  to  obtain  flow 
field  velocity  measurements.  A  high-frequency  tracking  shaft  angle 
encoder  was  used  to  divide  the  circumference  of  the  rotor  into  3000 
individual  bins.  The  circumferential  position  of  the  rotor  was  recorded 
with  each  velocity  measurement.  Each  of  the  15  main  blade  passages 
on  the  rotor  was  divided  into  200  bins  (N^p  =  200)  giving  a  pitchwise 
resolution  of  0.5  percent  of  main  blade  pitch.  Fiber  optic  probes  were 
used  to  track  blade  passing  so  that  the  shaft  angle  encoder  frequency 
could  be  reset  after  each  revolution.  The  encoder  count  at  the  end  of 
a  revolution  was  compared  to  3000.  If  the  difference  was  greater  than 
20  counts,  all  data  from  the  revolution  were  discarded  immediately. 
In  the  retained  data,  77  percent  of  the  revolutions  had  a  difference  of 
zero  counts  and  98  percent  were  within  5  counts  of  3000.  The  basic 
operating  principles  of  this  type  of  system  are  described  by  Strazisar 
(1986). 

Figure  4  shows  the  measurement  locations  in  the  impeller  and 
vaneless  diffuser.  Measurements  from  the  numbered  locations  will  be 
presented  herein.  The  spanwise  survey  lines  are  "quasi-orthogonal’'  or 
nearly  normal  to  the  hub  and  shroud  contours.  Control  of  the  viewing 
angle,  in  the  plane  of  Figure  4,  was  accomplished  by  mounting  the 
final  beam  projection  mirror  on  a  goniometric  cradle. 
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Impeller  Leading  Edge 


Figure  4:  Laser  Survey  Locations  in  4:1  Centrifugal 
Compressor 

Measurements  were  obtained  through  a  series  of  windows  located 
in  interchangeable  plugs  in  the  compressor  shroud.  Flat  windows  were 
used  to  eliminate  the  lensing  effects  which  arise  from  curved  windows. 
The  window  sizes  were  kept  small  to  minimize  the  local  clearance 
which  arises  from  the  mismatch  between  the  curved  shroud  and  the  flat 
windows. 

The  facility  was  seeded  by  atomizing  a  liquid  dispersion  of 
Alumina  (Aluminum-oxide)  powder  at  the  entrance  to  the  plenum  tank, 
far  upstream  of  the  measurement  locations.  The  powder  was 
suspended  in  100  percent  pure  ethanol.  The  electric  charge  on  the 
surface  of  the  Alumina  particles  in  solution  was  controlled  to  prevent 
particle  agglomeration.  Surface  charge  control  was  achieved  by 
adjusting  the  pH  of  the  ethanol-alumina  solution  using  small  amounts 
of  HCl  as  described  by  Wemet  and  Wemet  (1994).  The  seed  solution 
was  injected  into  the  plenum  tank  through  air  actuated  spray  nozzles 
which  produced  an  aerosol.  The  ethanol  component  evaporated 
completely  while  in  the  plenum  tank  leaving  only  Alumina  particles  to 
enter  the  compressor. 

The  Alumina  powder  used  in  this  investigation  had  a  specific 
gravity  of  3.96  and  a  mean  particle  size  of  0.7  ±  0.2  pm.  To  gage  the 
effect  of  the  relatively  high  specific  gravity  of  the  Alumina  powder  on 
velocity  measurement  accuracy,  a  back-to-back  test  was  done  in  a 
transonic  fan  using  both  Alumina  and  PSL  (poly-styrene-latex) 
particles.  The  PSL  used  in  the  comparison  had  a  mean  particle  size  of 
0.5  pm  and  a  specific  gravity  of  1.05.  The  results  obtained  with  each 
seed  material  were  nearly  identical  with  the  only  difference  occurring 
across  the  shock.  Since  the  flow  fields  in  the  current  investigation 
were  subsonic  throughout,  there  were  no  velocity  gradients  severe 
enough  to  effect  the  accuracy  of  velocities  measured  with  the  Alumina 
seed  material. 

Velocity  Determination 

The  quasi-meridional  survey  planes  containing  the  data  presented 
herein  are  planes  that  are  normal  to  the  plane  of  Figure  4  and  lie  along 


the  spanwise  survey  direction  shown  at  each  chordwise  location.  The 
laser  anemometer  acquired  velocity  measurements  in  a  measurement 
plane  that  was  orthogonal  to  both  the  spanwise  survey  direction  and 
the  plane  of  Figure  4.  A  rotating  beam  separation  prism  was  used  to 
rotate  the  laser  anemometer  fringe  pattern  in  the  measurement  plane. 
This  permitted  successive  measurements  of  several  components  of  the 
flow  vector  in  that  plane.  Components  of  the  absolute  velocity 
measurements  were  recorded  at  angles  of  0  and  +  15  degrees  away 
from  the  absolute  flow  angle  that  was  estimated  for  flows  in  the 
measurement  plane.  Spanwise  velocities  could  not  be  measured. 

For  each  probe  volume  orientation,  100,000  velocity 
measurements  were  typically  collected.  Data  collected  over  the  entire 
rotor  circumference  were  first  ensemble  averaged  into  a  single  blade 
passage.  The  resultant  flow  vector  (V^^s  p)  in  the  measurement 
plane  was  then  determined  for  each  measurement  bin  across  the  pitch 
by  fitting  the  ensemble-averaged  component  velocities  measured  at  the 
three  orientation  angles  to  a  cosine  function.  A  non-linear  least 
squares  technique  was  used  to  perform  the  fit. 

To  calculate  quasi-meridional  velocities  (VqJ  the  absolute 
velocity  at  each  measurement  location  was  multiplied  by  the  cosine 
of  the  absolute  flow  angle  measured  at  the  same  point.  The  pitchwise 
profiles  of  quasi-meridional  velocities  were  then  smoothed  using  a 
five  point  box-car  averaging  technique.  Equation  (1)  shows  the 
calculation  of  and  the  smoothing  technique. 
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Measurement  Uncertainty 

The  component  velocities,  measurement  angles,  and  component 
velocity  error  estimates  were  used  in  the  least  squares  fitting 
procedure  to  determine  the  velocity  vector  magnitude  and  flow  angle 
at  each  pitchwise  location.  The  least  squares  procedure  computes  the 
uncertainty  estimates  of  the  fitted  parameters  (velocity  and  flow  angle) 
based  on  the  sum  of  the  squared  residuals  of  the  fit  to  the  data,  the 
number  of  degrees  of  freedom  in  the  system  and  the  value  of  the 
parameter  element  in  the  covariant  matrix  (Wolberg,  1967). 

The  number  of  component  velocity  measurements  obtained  at 
each  fringe  pattern  orientation  were  not  distributed  evenly  across  a 
passage.  More  measurements  were  obtained  at  the  center  of  the 
passage  than  near  the  blade  surfaces.  The  uncertainty  estimate 
returned  by  the  least  squares  procedure  contains  the  effect  of  flow 
unsteadiness  and  blade-to-blade  geometry  variations  as  well  as  the 
effect  of  random  errors  in  measurement  technique. 

Plots  of  the  uncertainty  estimate  from  the  least  squares  procedure 
were  examined  at  each  measurement  plane.  The  uncertainty  in 
resultant  velocities  ranged  from  less  than  2%  in  the  center  of  a 
passage  to  30%  or  more  near  the  tip  and  blade  surfaces.  However, 
the  regions  of  high  uncertainty  coincide  with  areas  where  one  would 
expect  the  flow  unsteadiness  to  also  be  high.  Since  the  uncertainty 
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estimate  is  low  in  regions  where  low  unsteadiness  is  expected,  it  is 
assumed  that  measurement  uncertainty  due  to  technique  is  less  than 
2%. 

Figure  5  shows  a  comparison  of  laser  measurements  to  the 
absolute  velocities  determined  from  total  pressure  measurements  taken 
in  the  vaneless  diffuser  at  a  radius  ratio  of  1.10  which  is  slightly 
upstream  of  rake  location  A  in  Figure  1.  The  total  pressure 
measurements  were  made  using  a  probe  which  was  fully  immersed  in 
the  diffuser  passage  for  the  entire  survey.  As  the  survey  progressed, 
the  probe  traveled  deeper  into  a  pocket  in  the  diffuser  hub.  With  this 
arrangement,  the  blockage  created  by  the  probe  was  constant  over  the 
entire  span  of  the  survey.  The  static  pressures  needed  for  the  velocity 
calculation  were  measured  on  the  vaneless  diffuser  hub  and  shroud  at 
the  radius  of  the  probe.  A  linear  distribution  of  static  pressure  across 
the  passage  was  assumed  in  making  the  velocity  calculations. 


Figure  5:  Comparison  of  veiocities  measured  in  the 
vaneless  diffuser  using  LDV  to  those  caicuiated  from  totai 
pressure  measurements. 

COMPUTATiONAL  ANALYSIS 

The  computational  results  for  the  flow  field  were  obtained  using 
the  Reynolds-averaged  Navier-Stokes  code  developed  by  Dawes 
(1988).  The  code  solves  the  equations  of  motion  in  cylindrical 
coordinates  in  integral  conservation  form  using  six-sided  control 
volumes  formed  by  a  simple  H-mesh.  The  basic  algorithm  as  described 
by  Dawes  (1988)  is  similar  to  a  two-step  Runge-Kutta  method  plus 
residual  smoothing.  A  combined  second  and  fourth  derivative  artificial 
viscosity  model  with  pressure  gradient  switching  is  used  to  eliminate 
spurious  "wiggles"  and  to  control  shock  capturing.  The  eddy  viscosity 
is  obtained  using  the  Baldwin-Lomax  (1970)  mixing  length  model.  Tip 
clearance  is  handled  by  gradually  decreasing  the  thickness  of  the  blade 
to  zero  at  the  blade  tip  and  enforcing  periodicity  in  the  tip  gap. 

Grid  generation  was  done  using  a  NASA-developed  code  which 
accounts  for  the  splitter  blades  and  provides  good  flexibility  in 
distribution  of  the  grid  points.  Spacing  of  the  grid  points  is  by 
geometric  expansion  with  a  maximum  expansion  value  limited  by  the 


grid  code  to  1.3.  Although  expansion  factors  were  generally  much  less 
than  1.3,  values  of  1.3  at  the  impeller  inlet  and  1.17  at  the  exit  were 
used  in  order  to  define  the  tip  clearance  distribution.  For  the  analysis 
presented  herein  107  grid  points  were  used  in  the  streamwise  direction 
(J  direction)  with  the  main  blade  leading  edge  at  J=26,  the  splitter 
blade  leading  edge  at  J=52,  and  the  trailing  edge  of  both  at  J=88. 
Grid  points  were  clustered  around  the  main  and  splitter  blade  leading 
edges.  The  trailing  edge  was  tapered  to  zero  thickness  over  two  grid 
points.  Grid  spacing  in  the  streamwise  direction  varied  from  about 
0.3%  of  chord  near  the  leading  edges  to  2.6%  between  the  main  and 
splitter  blade  leading  edges,  to  nearly  5%  at  the  splitter  mid-chord. 
In  the  pitchwise  direction  45  grid  points  were  used  with  an  equal 
number  for  each  splitter  channel.  In  the  spanwise  direction  47  points 
were  used  with  4  cell  volumes  in  the  clearance  gap. 

The  grid  code  was  used  to  incorporate  the  measured  distribution 
of  tip  clearance,  which  varied  from  the  impeller  inlet  to  exit,  into  the 
grid.  The  Dawes  code  assumes  a  tip  clearance  that  is  a  constant 
percent  of  blade  span  from  inlet  to  outlet,  therefore,  it  was  modified 
to  read  the  grid  that  was  generated  to  match  the  measured  tip 
clearance  distribution.  The  full  tip  clearance  as  noted  under 
"Compressor  Performance"  was  used  in  the  analysis. 

The  Dawes  code  automatically  uses  the  wall-function 
approximation  for  the  wall  shear  stress  when  the  cell-center  nearest  a 
wall  lies  farther  from  the  surface  than  a  Y+  of  10.  For  the  solution 
presented  herein,  the  average  Y+  values  for  the  near-wall  cell-centers 
of  the  grid  were  19  for  the  hub,  32  for  the  shroud,  and  approximately 
60-75  for  all  blade  surfaces.  The  wall-function  was  therefore  used  at 
nearly  all  grid  points. 

A  number  of  grid  parametric  variations  were  done  to  investigate 
the  grid  sensitivity  of  the  solution.  The  effect  of  grid  density  was 
checked  by  increasing  the  pitchwise  number  of  grid  points  to  69  while 
maintaining  the  near-wall  spacing,  which  reduced  the  grid  spacing  by 
50%  at  mid-pitch.  The  number  of  grid  points  in  the  spanwise 
direction  was  increased  to  77  while  maintaining  the  near- wall  spacing, 
which  reduced  the  spanwise  grid  spacing  by  50%.  These  changes  in 
the  grid  had  a  relatively  small  effect  on  the  solutions. 

The  sensitivity  of  the  solution  to  the  treatment  of  the  tip 
clearance  was  checked  by  generating  a  grid  with  0.6  times  the 
physical  clearance  and  another  grid  with  the  inlet  tip  clearance 
doubled.  The  number  of  cells  in  the  tip  clearance  was  also  doubled 
from  4  to  8.  These  changes  had  a  small  effect  on  the  solution  and  no 
effect  on  the  conclusions  reached  in  this  work.  The  amount  of 
clearance  flow  that  crossed  the  blade  tips  was  calculated  to  be  7.8% 
of  the  inlet  mass  flow  for  the  main  blade  and  5.8%  for  the  splitter 
blade  using  the  solution  generated  with  the  baseline  grid. 

RESULTS  AND  DISCUSSION 
Inlet  Region 

The  operating  point  selected  for  all  laser  anemometer  surveys 
and  computations  reported  here  was  the  design  mass  flow  rate  of  4.54 
kg/s  (10  Ibm/s)  with  the  previously  described  tip  clearance 
distribution.  This  point  is  highlighted  in  Figure  3. 
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The  normalized  through-flow  velocity  distribution,  Vq^/Utip, 
measured  at  1%  of  main  blade  chord  is  shown  in  a  carpet  plot  format 
in  Figure  6.  Since  velocity  magnitudes  are  difficult  to  determine  from 
a  three  dimensional  plot,  a  table  is  included  to  give  a  snapshot  of  the 
velocity  magnitudes.  Spans  are  listed  from  tip  to  hub,  where  100 
percent  span  is  located  at  the  blade  tip.  At  each  span,  the  normalized 
through-flow  velocity  is  reported  at  25,  50  and  75  percent  blade  gap 
across  the  passage.  The  velocity  profiles  at  each  spanwise  location 
contain  a  small  at  points  where  the  velocity  magnitudes  are 
reported.  The  velocity  distribution  shown  in  Figure  6  increases  toward 
the  suction  surface  in  a  linear  fashion  and  is  similar  to  inducer  velocity 
profiles  that  were  reported  by  Eckardt  (1976),  Krain  (1987),  and 
Hathaway,  et  al.  (1993). 


Figure  6:  Distribution  of  through-flow  velocities  measured 
on  a  plane  located  at  1%  of  main  blade  chord 


Figures  7-9  show  carpet  plots  of  data  taken  at  5%,  10%,  and 
29%  of  main  blade  chord,  respectively.  All  plots  show  the  typical 
increase  in  velocity  from  pressure  surface  to  suction  surface  with  some 
indication  of  low  momentum  flow  near  the  shroud.  This  low 
momentum  region  is  first  detected  near  the  suction  side  of  the  channel 
at  5%  chord  and  100%  span.  At  10%  chord  the  velocity  in  this  region 
has  decreased  considerably  relative  to  the  velocity  at  95%  span.  The 
plane  at  29%  chord  does  not  show  any  evidence  of  the  low  momentum 
region.  Data  was  not  acquired  beyond  95%  span  on  this  plane,  and  the 
low  momentum  region  near  the  tip  may  well  lie  outboard  of  95%  span. 
The  flow  field  at  29%  chord  does,  however,  show  evidence  of  the 
potential  field  of  the  splitter  leading  edge  located  at  30%  chord. 

Data  acquired  downstream  of  the  splitter  leading  edge  will  be 
presented  in  the  form  of  cross-passage  contour  plots  of  Vqj/U^p  since 
contour  plots  illustrate  the  migration  of  low  velocity  regions  better  than 
carpet  plots.  In  the  discussion  of  these  results,  the  passage  adjacent  to 
the  main  blade  pressure  surface  will  be  referred  to  as  the  "pressure 
side"  passage.  The  "suction  side  passage"  will  refer  to  the  passage 
adjacent  to  the  main  blade  suction  surface. 


Figure  7:  Distribution  of  through-flow  velocities  measured 
on  a  plane  located  at  5%  of  main  blade  chord 
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Figure  8:  Distribution  of  through-flow  velocities  measured 
on  a  plane  located  at  10%  of  main  blade  chord 


Carpet  and  contour  plots  of  the  through-flow  velocity  distribution 
at  the  splitter  leading  edge  plane  at  30%  of  main  blade  chord  are 
shown  in  Figures  10  and  11.  The  carpet  plot  shows  the  obvious 
presence  of  the  splitter  but  also  indicates  a  large  region  of  low 
momentum  flow  near  the  tip  on  either  side  of  the  splitter  leading  edge. 
This  is  shown  more  clearly  in  the  contour  plot,  where  the  loading  on 
the  splitter  leading  edge  can  be  seen  at  aU  spans.  Low  momentum 
regions  occur  around  the  splitter  leading  edge  and  near  the  main  blade 
pressure  surface.  The  latter  results  from  deceleration  of  the  flow  in 
the  streamwise  direction  along  the  pressure  surface  due  to  the  main 
blade  loading. 
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Figure  9:  Distribution  of  through-flow  velocities 
measured  on  a  plane  located  at  29%  of  main  blade 
chord 


Figure  10;  Distribution  of  through-flow  velocities  on  a 
plane  located  near  the  splitter-blade  leading  edge  (30% 
of  main-blade  chord) 


A  clearer  picture  of  the  flow  development  upstream  of  the  splitter 
leading  edge  can  be  obtained  using  the  simulation  generated  with  the 
Dawes  code.  Figure  11a  shows  the  predicted  distribution 
on  cross-sectional  planes  at  5,  10  and  30%  chord,  properly  oriented  in 
the  impeller  channel.  The  calculations  at  these  three  planes  show 
qualitatively  the  same  information  as  did  the  measured  data.  The 
predicted  momentum  deficit  near  the  blade  tip  is  somewhat  more 


Figure  11:  Through-flow  velocity  contours  on  a  plane 
located  near  the  splitter-blade  leading  edge  (  30%  of 
main  blade  chord) 


shallow  and  broader  at  5%  chord  than  at  10%  chord  (note  the 
velocities  at  50%  and  75%  pitch  at  100%  span  in  figures  7  and  8). 
The  tip  clearance  flow  impacts  the  splitter  leading  edge  as  it  moves 
across  the  main  blade  channel  from  the  suction  side  toward  the 
pressure  side.  Relative  velocity  vectors,  not  shown  here,  were 
generated  on  the  grid  surfaces  just  above  the  tip  of  the  blade  and  were 
color  coded  to  match  the  quasi-meridional  velocity  shown  on  the 
cross-sections.  These  vectors  clearly  showed  the  movement  of  the 
clearance  flow  across  the  channel.  Flow  tracers  released  in  the  tip  gap 
near  the  suction  side  of  the  passage  at  several  percent  chord  locations 
are  shown  in  Figure  11a  to  indicate  this  general  movement.  The  red 
trace  which  was  released  on  the  suction  surface  near  the  leading  edge 
migrates  radially  outward  to  the  blade  tip  where  it  is  entrained  by  the 
clearance  flow  and  pushed  toward  the  splitter. 

The  secondary  flow  in  the  inlet  is,  as  might  be  expected  for  a 
subsonic  impeller,  quite  similar  to  that  reported  by  Hathaway  et  al. 
(1993)  for  a  large  centrifugal  compressor  in  which  quite  detailed 
three-dimensional  velocities  were  measured  with  a  laser  anemometer 
system.  The  particle  tracer  studies  of  the  CFD  simulation,  not  shown 
here,  indicate  that  fluid  near  the  blade  and  hub  solid  surfaces  is  driven 
toward  the  shroud  by  the  hub-to-shroud  pressure  gradient  established 
by  the  main  flow.  Once  at  the  shroud  this  low  momentum  flow  is 
entrained  by  the  flow  leaking  over  the  tip  of  the  blade  and  carried 
toward  the  pressure  side  of  the  channel.  This  was  clearly  shown  in 
the  work  of  Hathaway  et  al.  with  both  the  measurements  and 
calculations.  The  only  difference  in  the  present  work  is  the  presence 
of  the  splitter  which  causes  the  deficit  flow  to  be  unequally  distributed 
in  the  two  channels  formed  by  the  splitter  blade. 
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Figure  11a:  Predicted  distribution  of  at  5,  10  and  30%  of  main  blade  chord 


Splitter  Entry  Region 

At  52%  chord,  Figure  12,  the  pitchwise  extent  of  the  velocity 
deficit  has  increased  in  both  passages.  The  magnitude  and  spanwise 
depth  of  the  deficit  region  have  increased  much  more  in  the  suction 
side  passage  where  the  depth  extends  down  to  80%  span.  The  center 
of  this  region  is  nearer  the  splitter  blade  pressure  surface.  Data  not 
shown  here,  that  were  collected  at  85%  and  90%  spans  on  a  plane  at 
47%  chord,  also  indicate  the  presence  of  a  well  developed  low  velocity 
region  near  the  splitter  blade  in  the  suction  side  passage.  The 
pitchwise  coverage  of  the  deficit  has  increased  in  the  pressure  side 
passage,  while  the  spanwise  depth  and  magnitude  of  the  velocity 
gradient  have  not  changed  very  much.  The  center  of  the  deficit  region 
appears  closer  to  the  center  of  the  passage  than  it  did  at  30%  chord,  as 
though  it  is  still  moving  toward  the  main-blade  pressure  surface. 


%m/^tip  OAS  0.20  0.24  0.28  0.32  0.36  0.40  0.44 


Figure  12:  Contours  of  through-flow  velocity  measured 
on  a  plane  located  at  52%  of  main  blade  chord 
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The  large  deficit  region  in  the  suction  side  passage  at  52%  chord 
results  when  the  splitter  "scrapes  up"  the  tip  clearance  and  near  wall 
flows  that  have  been  moving  across  the  main  blade  channel  from  the 
suction  surface  toward  the  pressure  surface.  This  process  is  more 
clearly  seen  in  Figure  12a  which  shows  flow  tracers  generated  from  the 
CFD  solution  superimposed  on  a  blade-to-blade  plot  of  the  quasi¬ 
meridional  velocity  measured  with  the  laser  anemometer  at  90%  span. 
The  black  tracers  released  in  the  main  blade  tip  clearance  gap  from  the 
leading  edge  to  74%  chord  are  swept  toward  the  pressure  side  of  the 
channel.  The  tracers  upstream  of  the  splitter  leading  edge  are 
continually  pushed  deeper  into  the  main  flow  by  the  leakage  flow 
occurring  downstream  and  are  subsequently  turned  in  the  streamwise 
direction  by  the  main  flow.  This  "early"  leakage  flow  is  "scraped"  by 
the  splitter  and  forms  the  low  momentum  region  on  the  splitter 
pressure  surface  as  shown  in  figures  10,  11  and  12.  As  the  fluid 
particles  in  this  leakage  flow  hit  the  splitter  they  are  turned  down  into 
the  main  flow  to  form  the  deficit  in  velocity  in  the  suction  side 
passage  at  52%  chord,  as  shown  in  Figure  12a. 


The  white  tracers  in  Figure  12a  are  released  along  the  span  near 
the  main  blade  suction  surface  at  the  splitter  leading  edge  plane. 
These  tracers  move  toward  the  shroud  and  are  entrained  by  the 
cross-channel  flow  near  the  shroud.  The  cyan  tracers  are  released  in 
the  tip  clearance  gap  downstream  of  the  splitter  leading  edge.  Some 
of  the  fluid  marked  by  both  the  cyan  and  white  tracers  is  forced  into 
the  deficit  region  in  the  suction  side  passage  and  some  combines  with 
the  leakage  from  the  main  blade  and  leaks  over  the  splitter  blade  and 
into  the  deficit  region  in  the  pressure  side  passage.  The  deficit  region 
in  the  pressure  side  passage  is  not  as  large  as  that  in  the  suction  side 
passage  because  the  leakage  flow  generated  from  the  main  blade 
upstream  of  the  splitter  leading  edge  is  primarily  contained  in  the 
suction  side  channel.  The  simulations  show  the  deficit  region  to  be 
closer  to  the  center  of  the  suction  side  channel  than  does  the  data. 
According  to  the  simulations,  52%  chord  is  just  downstream  of  where 
the  low  momentum  fluid  impacts  on  the  splitter  blade  pressure 
surface  and  begins  to  move  back  toward  the  suction  surface  of  the 
main  blade. 


Figure  12a:  CFD  tracers  superimposed  on  measured  through-flow  velocity  at  90%  span 
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Exit  Region 

The  measured  flow  field  downstream  of  70%  chord  is  presented 
in  Figures  13-15,  17  and  18.  The  predicted  flow  field  at  70,  84  and  96 
percent  chord  is  shown  in  Figure  16. 

By  70%  chord  (Figure  13)  the  shroud  deficit  region  in  the  suction 
side  passage  has  moved  closer  to  the  suction  surface  of  the  main  blade 
and  a  deficit  has  also  developed  in  the  pressure  side  passage  near  the 
splitter  blade  suction  surface.  The  low  momentum  regions  near  the 
pressure  side  of  each  passage  are  due  to  the  loading  on  the  main  and 
splitter  blades  and  are  comparable  in  size.  By  this  location,  the 
calculations  have  begun  to  deviate  more  from  the  data  but  the  general 
flow  pattern  is  captured  (i.e.  there  is  more  low  momentum  fluid 
accumulated  near  the  shroud  for  the  suction  side  passage  than  for  the 
pressure  side  passage  in  both  the  data  and  the  calculations).  The 
calculation,  however,  indicates  that  a  low  momentum  “tail"  is 
beginning  to  form  and  is  being  pushed  toward  mid-span  in  both 
passages  as  shown  in  Figure  16.  That  feature  is  not  seen  in  the  data. 
The  velocity  deficit  for  the  calculations  is  greater  than  that  of  the  data. 


released  in  the  calculations  in  the  center  of  the  "tail"  of  low  velocity 
fluid  in  the  pressure  side  passage,  show  that  this  fluid  came  mainly 
from  boundary  layer  fluid  on  the  main  blade  pressure  surface  but  with 
some  contribution  from  the  migration  of  boundary  layer  fluid  on  the 
splitter  blade  suction  surface  and  tip  clearance  flow  over  the  splitter 
blade. 


Figure  14:  Contours  of  through-flow  velocity  measured 
on  a  plane  located  at  84%  of  main-blade  chord 
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Figure  13:  Contours  of  through-flow  velocity  measured  on 
a  plane  located  at  70%  of  main-blade  chord 


At  84%  chord,  Figure  14,  the  low  velocity  region  is  now 
accumulating  near  the  suction  surface  in  both  passages  and  is 
increasing  in  depth.  In  the  pressure  side  passage,  a  low  velocity  core 
has  moved  downward  in  spanwise  position  near  the  splitter  blade.  This 
core  now  contains  the  lowest  through-flow  velocity  fluid  found  in 
either  passage.  The  velocity  distribution  found  at  84%  chord  is  the 
opposite  of  that  found  farther  upstream  at  52%  chord  (Figure  12) 
where  the  suction  side  passage  contained  the  lowest  velocity  flow.  In 
the  suction  side  passage  at  84%  chord  the  low  velocity  region  near  the 
main  blade  suction  surface  has  the  same  general  shape  as  the  region 
in  the  pressure  side  passage  but  does  not  contain  a  well  developed 
core  of  very  low  velocity  fluid.  The  low  velocity  region  near  the 
splitter  blade  in  the  pressure  side  channel  is  quite  similar  to  the  low 
velocity  region  in  the  pressure  side  passage  for  the  calculations.  The 
calculations  (Figure  16)  show  very  low  velocity  fluid  near  the  shroud 
with  the  "tail"  mentioned  at  70%  chord  more  pronounced  and  closer 
to  50%  span.  Although,  there  does  not  appear  to  be  a  "tail"  in  the 
suction  side  passage  for  the  data,  with  suitable  high-lighting  of  the 
contours,  a  faint  tail  is  distinguishable.  It  is  not,  however,  as 
pronounced  as  the  one  produced  by  the  calculations.  Flow  tracers 


At  96%  cord.  Figure  15,  the  velocity  deficit  remains 
well-developed  in  the  pressure  side  passage.  A  large  deficit  region 
exists  near  the  shroud  from  about  50%  pitch  to  the  splitter  suction 
surface  and  at  about  50%  span  near  the  same  surface.  The  deficit  in 
the  pressure  side  passage  looks  somewhat  similar  to  that  produced  by 
the  calculations  but  the  velocity  magnitudes  are  different.  The  flow 
pattern  measured  in  the  suction  side  channel  is  similar  to  that  in  the 
pressure  side  but  the  velocity  gradients  are  far  less. 


Vqm/^tip  0-20  0-23  0.26  0.29 


Figure  15:  Contours  of  through-flow  velocity  measured  on 
a  plane  located  at  96%  of  main  blade  chord 

Based  upon  the  calculations  and  their  agreement  with  the  data 
we  see  a  similar  type  of  flow  in  the  pressure  and  suction  side 
passages.  Low  momentum  tip  clearance  flows  and  blade  viscous 
layers  migrate  in  the  channel  and  accumulate  near  the  shroud-suction 
surface  in  both  passages  which  is  consistent  with  observations  by 
other  authors.  It  is  not  clear,  however,  why  the  deficit  region  in  the 
pressure  side  passage  is  more  pronounced  than  that  in  the  suction  side 
passage.  Integration  of  the  calculations  to  determine  the  mass  flow 
split  in  each  splitter  channel  revealed  that  the  flow  in  the  pressure  side 
passage  was  consistently  0.5  to  1%  greater  than  that  in  the  suction 
side  passage. 
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Figure  16:  CFD  prediction  of  through-flow  velocity  at  70,  84  and  96  percent  of  main  blade  chord 


At  99%  chord,  Figure  17,  the  measured  through-flow  velocity 
contours  have  taken  on  a  more  stratified  appearance  at  pitchwise 
locations  that  are  away  from  the  suction  surface  in  each  passage. 
Pockets  of  low  velocity  flow  remain  centered  about  a  lower  span  in  the 
region  adjacent  to  each  suction  surface. 
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Figure  17:  Contours  of  through-flow  velocity  measured 
on  a  plane  located  at  99%  of  main  blade  chord 

Except  for  the  pockets  of  low  velocity  fluid  from  Figure  17,  the 
profile  across  each  passage  near  the  discharge  is  somewhat  linear 
across  the  pitch  and  does  not  resemble  a  traditional  jet-wake  profile. 
Figure  18  more  clearly  shows  the  through-flow  velocity  profile  at  99% 
chord.  The  low  velocity  pockets  of  Figure  17  can  not  be  seen  as  they 
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Figure  18:  Distribution  of  through-flow  velocities  on  a 
plane  located  at  99%  of  main  blade  chord 
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are  behind  the  profiles  shown  for  the  outermost  spans.  In  the  pressure 
side  passage,  the  profile  is  nearly  linear  across  the  pitch.  The  profile 
of  the  suction  side  passage  is  similar  to  the  discharge  profile  reported 
by  Krain  (1987)  for  an  impeller  with  30  degrees  of  backsweep.  The 
discharge  profile  of  that  impeller  showed  the  through-flow  velocity 
increasing  toward  the  hub  and  less  so  toward  the  pressure  surface  with 
a  slight  depression  near  midpitch. 

Wake  Decay  in  the  Vaneless  Diffuser 

Hathaway  et  al.  (1993)  reported  wake  decay  from  their  large 
centrifugal  compressor  tests  which  indicated  that  the  wake  from  the 
blades  decays  very  rapidly  such  that  at  radius  ratios  of  1 .04  and  1 .06 
the  pitchwise  variation  in  absolute  flow  angle  is  below  10  degrees.  In 
addition  they  measured  the  decay  of  the  wake  near  the  shroud  where 
the  tip  clearance  flow  formed  a  second  region  of  low  radial  momentum 
fluid.  The  blade  wake  fluid  near  the  shroud  also  decayed  very  rapidly 
near  the  trailing  edge;  however,  the  "wake"  associated  with  clearance 
flow  decayed  less  rapidly  as  a  result  of  its  much  lower  velocity 
gradient  and,  thus,  shear  stress  compared  to  that  in  the  blade  wakes. 

In  the  present  investigation  spanwise  surveys  were  taken  at  radius 
ratios  of  1.01,  1.02,  1.04,  1.06,  1.08,  and  1.10.  The  angle  variation  is 
shown  in  Figure  19.  The  wakes  have  been  shifted  in  the  pitchwise 
direction  to  align  all  of  the  main  blade  wakes  in  the  plot.  Each  trace 
is  referenced  to  the  minimum  flow  angle,  radius  ratio  and 

percent  span  in  order  to  show  the  pitchwise  variation  to  which  a 
diffuser  vane  might  be  subjected  if  one  were  present.  The  angle 
variation  decreases  rapidly  to  10  degrees  or  less  by  a  radius  ratio  of 
about  1 .04  which  is  consistent  with  the  value  reported  by  Hathaway 
for  comparable  radius  ratios.  In  the  present  study  the  tip  clearance 
"wake"  is  adjacent  to  the  blade  wake  suction  surface  and  is 
distinguished  from  the  blade  wake  at  and  above  85%  span  by  a  distinct 
"plateau"  in  absolute  flow  angle  over  about  15%  of  main  blade  pitch 
at  the  100%  percent  span  location.  This  is  due  to  the  clearance  flow 
accumulating  near  the  blade  surfaces  at  the  impeller  exit  as  indicated 
in  Figure  17. 


%  Pitch 

Figure  19:  Pitchwise  absolute  flow-angle  difference  near 
100  percent  of  diffuser  span 


CONCLUSIONS 

1)  The  formation  of  a  through-flow  velocity  deficit  may  begin  farther 
upstream  than  has  been  previously  suggested  for  high  speed  impellers. 
It  begins  as  the  inducer  suction  surface  boundary  layer  migrates 
toward  the  blade  tip  and  becomes  entrained  in  the  tip  clearance  flow. 
The  deficit  then  moves  toward  the  pressure  surface. 

2)  The  majority  of  the  initial  deficit  is  trapped  and  "scraped  up"  by 
the  splitter  blade  which  causes  it  to  turn  down  and  roll-up  in  the 
suction  side  passage,  thereby  creating  a  larger  velocity  deficit  in  that 
passage. 

3)  Some  flow  leaks  over  the  splitter  blade  tip  in  a  manner  that  is 
similar  to  the  inducer  leakage  flow.  This  flow  joins  with  the  portion 
of  the  upstream  deficit  that  entered  the  pressure  side  passage. 

4)  The  velocity  deficits  in  both  passages  are  eventually  driven  back 
toward  the  suction  surface  by  the  cross-passage  pressure  gradients. 

5)  The  discharge  velocity  profile  of  this  impeller  did  not  have  the 
characteristics  of  a  jet  wake.  The  profile  was  similar  to  discharge 
profiles  reported  by  others  who  have  tested  backswept  impellers. 

6)  The  rate  of  decay  of  the  overall  wake  region  is  similar  to  that 
which  was  reported  by  Hathaway  et  al.  (1993)  in  a  large  low-speed 
impeller. 
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